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ABSTRACT 

We investigate the formation and evolution of interstellar dust-grain ices under dark-cloud con- 
ditions, with a particular emphasis on CO2. We use a three-phase model (gas/surface/mantle) to 
simulate the coupled gas-grain chemistry, allowing the distinction of the chemically-active surface 
from the ice layers preserved in the mantle beneath. The model includes a treatment of the compe- 
tition between barrier-mediated surface reactions and thermal-hopping processes. The results show 
excellent agreement with the observed behavior of CO2, CO and water ice in the interstellar medium. 
The reaction of the OH radical with CO is found to be efficient enough to account for CO2 ice pro- 
duction in dark clouds. At low visual extinctions, with dust temperatures >^12 K, CO2 is formed by 
direct diffusion and reaction of CO with OH; we associate the resultant C02-rich ice with the obser- 
vational polar C0 2 signature. CH 4 ice is well correlated with this component. At higher extinctions, 
with lower dust temperatures, CO is relatively immobile and thus abundant; however, the reaction 
of H and O atop a CO molecule allows OH and CO to meet rapidly enough to produce a CO:C02 
ratio in the range ~2— 4, which we associate with apolar signatures. We suggest that the observational 
apolar CO2/CO ice signatures in dark clouds result from a strongly segregated CO:H20 ice, in which 
CO2 resides almost exclusively within the CO component. Observed visual-extinction thresholds for 
CO2, CO and H 2 are well reproduced by depth-dependent models. Methanol formation is found to 
be strongly sensitive to dynamical timescales and dust temperatures. 

Subject headings: Astrochcmistry, ISM: abundances, ISM: clouds, ISM: molecules, molecular processes 



1. INTRODUCTION 

Interstellar ices are detected in both quiescent and star- 
forming regions in the ISM, through infrared absorption 
line spectroscopy. These ices exist as mantles around 
a dust-grain core, and are believed to form by the ac- 
cretion of gas-phase atomic and molecular material onto 
the grain surfaces, whereupon further chemical process- 
ing occurs. Observational data provide a number of con- 
straints both on the chemical compositions of these ices 
and on their physical structure. Water ice constitutes 
the primary component of interstellar ice, but significant 
quantities of other molecules are also present, including 
CO, C0 2 and methanol (CH 3 OH) (Chiar et al. 1995, 
Whittet et al. 2001, 2007), with abundances with respect 
to water ranging into the tens of percent, dependent on 
the astronomical source (Gibb et al. 2010). Comparison 
of observational data with laboratory spectra may also 
indicate the nature of the ice mixture in which the origi- 
nators of the line profiles reside, i.e. whether dominated 
by polar water ice or by a non-polar mixture in which 
H2O is a minor component. 

Observations indicate that water is the first ice con- 
stituent to attain abundances sufficient to be detected in 
the ISM; the most recently obtained threshold for H2O- 
ice detection (using the 3 /xm line) lies at a visual extinc- 
tion of 3.2 ± 0.1 (Whittet et al. 2001), indicating that 
the formation/deposition of this molecule occurs under 
conditions in which the interstellar radiation field (ISRF) 
is attenuated by a visual extinction of at least 1.6 (repre- 
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senting the extinction from the edge to the center of the 
cloud). The observed (edge-to-edge) visual extinction 
threshold for CO2 of 4.3 ± 1.0 is close to that of water ice, 
while CO has a significantly greater threshold of 6.7± 1.6 
(Whittet et al. 2007). Gas-phase formation of CO2 is 
far less efficient than that of CO, thus grain-surface CO 
(which is itself primarily formed in, and accreted from, 
the gas phase) is believed to be the progenitor of CO2 
ice. The prevalence of CO2 ice, and lack of CO ice, at 
low extinctions requires the CO-to-C02 conversion to be 
highly efficient. However, a significant drop in the effi- 
ciency at higher extinctions is then required to account 
for the detection of CO ice itself. It is not clear whether 
the changes in chemical behavior are caused directly by 
the differing visual extinctions, or whether they are re- 
lated to the gas dynamics of the clouds, or to some other 
time-dependent effect. 

The problem is further complicated by the fact that 
there is no clear, efficient, chemical pathway that can 
account for the majority of CO2 ice formation. Three 
processes are typically cited as the most likely: 

HCO + O^C0 2 +H (1) 

CO + OH^C0 2 +H (2) 

CO + O -> C0 2 (3) 

Reaction (1) is barrierless, reaction (2) is typically as- 
signed a small activation energy barrier (80 K, Ruffle 
& Herbst, 2001), and reaction (3) is thought to require 
a more substantial activation energy of 290 - 1000 K 
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(Roser et al. 2001; d'Hendecourt et al. 1985). The re- 
actants O and OH are thought to be prevalent on the 
grains, as they are involved in the formation of water 
ice. However, the HCO radical required for reaction (1) 
is formed by the hydrogcnation of CO on the grains, 
and this process has an activation energy barrier of at 
least 500 K (Woon, 2002). Chemical models of inter- 
stellar clouds typically find rather that it is reaction (3) 
that shows the most promise for CO2 formation; in con- 
trast to HCO, CO is expected to attain at least moderate 
abundances on the grain surface, even if CO2 formation 
is reasonably efficient. Reaction (2) is usually found to 
be hampered by the relative immobility of CO and OH 
at low temperatures; hydrogenation of these species is 
far more competitive. However, Oba et al. (2010) re- 
cently conducted surface experiments to isolate reaction 
(2), allowing CO to react with the non-energetic prod- 
ucts of H 2 photolysis, at temperatures as low as 10 K. 
The process was found to produce detectable quantities 
of C0 2 . 

Tielens & Hagen (1982) suggested that the formation 
of CO2 on grains requires reaction (3) to proceed with 
high efficiency, and later studies such as Ruffle & Hcrbst 
(2001) and Hassel et al. (2010) have found that only 
models with severely reduced activation energy barriers 
for reaction (3) are capable of reproducing observed CO2 
ice abundances; atomic oxygen is otherwise too quickly 
hydrogenated before reaction (3) can occur. Even such 
low-barrier models cannot apparently reproduce the ex- 
tinction dependence of CO and C0 2 formation. More 
recently, Charnley & Rodgers (2009) demonstrated, us- 
ing a stochastic three-phase chemical model, that reac- 
tion (3) could not permit the co-existence of polar and 
apolar C0 2 ices on a single grain; these authors instead 
proposed reaction (1) as the primary formation mecha- 
nism. 

Goumans & Andersson (2010) recently calculated the 
temperature-dependent rate of reaction (3) in the gas 
phase. They found the rate to be yet lower than is 
typically assumed in models. However, they also sug- 
gested that CO and O may form a loosely-bound O...CO 
complex on the grain surface; the hydrogcnation of the 
oxygen atom in the complex would immediately lead to 
the formation of a highly-excited OH. ..CO system that 
could easily overcome the moderate barrier to reaction 
(2). This mechanism has yet to be tested in computa- 
tional models. 

The formation of simple dust-grain ices now appears 
to be crucial to the formation of much more com- 
plex molecules observed in the gas phase toward star- 
forming regions. Ubiquitous molecules such as methanol 
(CH3OH) and methyl formate (HCOOCH3) exhibit in- 
terstellar abundances much greater than may be ac- 
counted for by gas-phase processes (e.g. Horn et al. 2004, 
Geppert et al. 2006, Garrod et al. 2007). Understanding 
the earliest phases of ice formation is therefore essential 
to our knowledge of the chemistry of later evolutionary 
stages of star formation. 

Here we employ a three-phase chemical model to dis- 
tinguish between the chemically-active surface layer and 
the deeper, inactive, mantle layers within the ice. The 
formation of the ice mantle allows the preservation of the 
chemical composition of earlier surface layers during the 
evolution of the cloud. By monitoring the deposition of 



mantle material over time, the model reveals the com- 
position of each layer within the mantle, allowing the 
specific ice mixtures within which important molecules 
reside to be assessed, and the relationship between ice 
mantle composition and earlier gas-phase chemistry to 
be determined. In this paper we analyse the modeled 
abundances of the main observationally-determined com- 
ponents of interstellar dust-grain ices, and examine their 
formation mechanisms and time-dependent behavior. 

In the following section we introduce the details of the 
chemical model and highlight changes made in the treat- 
ment of chemical processes. In section 3 we present the 
results of the model, applied under a variety of physical 
conditions. A discussion of the results follows in section 
4. The conclusions are summarized in section 5. 

2. METHOD 

We use a gas-grain chemical model to investigate the 
time-dependent evolution of chemical abundances in a 
cold dark interstellar cloud. The basic model is that used 
by Garrod et al. (2008), and includes the grain-surface 
and gas-phase reactions used by those authors, as well 
those added to the network by Bcllochc et al. (2009). 
The gas-phase reaction rates not associated with the ad- 
dition of new complex molecules by the aforementioned 
authors derive from the osu.2005 chemical network. The 
network includes approximately 5700 gas-phase reac- 
tions, and 2900 surface/mantle processes including evap- 
oration, accretion, and surface-mantle transport. 

We employ the grain-surface binding/desorption en- 
ergies (-Edos) an d diffusion barriers (i?dif) of Garrod & 
Herbst (2006) and those values added by Garrod et al. 
(2008), all of which are representative of a surface com- 
posed primarily of amorphous water-ice. Table 1 shows 
a selection of values for the most important surface re- 
actants. The diffusion barrier of each surface species, i, 
is uniformly fixed to a fraction of the associated bind- 
ing energy, initially set to E^f(i) — 0.5 E dcs (i). Values 
ranging from 0.3 to 0.8 E des (i) have been adopted in the 
past (Tielens & Hagen 1982; Katz et al. 1999; Ruffle & 
Herbst 2000; Cuppen et al. 2009); however, within this 
range, the values are generally very poorly constrained. 

The basic Ed cs and i?dif values are further adjusted 
according to the composition of the ice surface - in par- 
ticular, the fraction of the ice composed of molecular 
hydrogen, H 2 . When assuming dust temperatures < 10 
K, under conditions where H 2 is abundant in the gas 
phase, gas-grain models can produce solutions in which 
the dust-grain ices are composed primarily of H 2 , and 
the problem is exascerbated by the use of a three-phase 
chemical model (see section 2.2). This behavior is un- 
physical, and is caused by the use of binding energies 
that are appropriate to a water-ice surface, rather than 
to a mixed H 2 :H 2 surface, which should demonstrate 
much weaker binding. Cuppen et al. (2009) estimate 
that the binding energy to an H 2 surface is around 10 
times weaker than that for CO. Since, as also explained 
by Cuppen et al., the total binding energy is an aggregate 
of contributions from neighboring surface constituents, 
we calculate effective binding energies and diffusion bar- 
riers according to the fractional coverage of the surface 
with H 2 , 0(H 2 ), giving: 

E dcs . cB = E dcs [1 - 0(H 2 )] + 0.1 E dcs 0(H 2 ) (4) 
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with the same relationship for £?dif,eff- So long as the 
modifying factor used in equation (4), here 0.1, is sig- 
nificantly less than unity then the precise value assumed 
should be of small consequence. Since molecular hydro- 
gen is the only chemical species likely to build up an 
appreciable grain-surface abundance that could alter the 
surface binding properties so strongly, as compared to the 
water-ice surface, we restrict these considerations only to 
the H 2 fraction. However, for completeness we allow the 
binding properties of all species to be altered according to 
this fraction, as per equation (4). This treatment renders 
.Edes an d £W as time-dependent quantities, making the 
differential equations somewhat stiffer than otherwise. A 
more comprehensive treatment of effective binding ener- 
gies may be attempted in future. The current method 
produces maximum H 2 -ice fractions of around 10% or 
less, under dark-cloud conditions in which molecular hy- 
drogen dominates in the gas phase. The variation of the 
binding energies tends to flatten out quickly as the H 2 
fraction stabilizes, according to the balance between ac- 
cretion and thermal evaporation; other processes have a 
minimal effect on surface H2 abundance. 

The grain-surface chemistry is treated using the 
modified-rate method of Garrod (2008), to account for 
stochastic effects where present. The method has been 
shown by Garrod et al. (2009) to produce an excellent 
match to exact Monte Carlo simulation methods, over a 
broad range of physical parameters. A canonical grain 
size of O.l^m is assumed, with 10 6 surface binding sites. 

The model includes gas-phase and grain-surface 
cosmic-ray photodissociation/photoionization processes 
(as per Garrod et al. 2008), assuming a cosmic-ray ion- 
ization rate of £ = 1.3 x 10~ 17 s _1 . The initial elemental 
abundances used in all models are shown in table 2, and 
are the same as those used by Garrod et al. (2008). 

We include the reactive desorption mechanism given 
by Garrod et al. (2007), whereby the formation of a 
molecule on the grain surface can result in its desorption 
into the gas phase, with a probability of jl %. We assign 
a value orrk = 0.01 for this process. We also include 
photodesorption for all surface species, using the rates 
and yields of Oberg et al. (2009a, b); yields of 10 -3 are 
assumed for species for which no data currently exist. 
For visual extinctions much larger than 1, reactive des- 
orption strongly dominates over direct photodesorption; 
although, in some cases, the reactive desorption process 
is preceded by the photodissociation of a surface molecule 
by the cosmic ray-induced UV field. 

2.1. Reaction-diffusion competition 

The efficiencies of activation-energy barrier-mediated 
reactions on the grain surfaces are crucial to this study. 
In general, under the classical rate-equation expression, 
the rate of a reaction A+B is expressed as: 

k AB = fact(AB) x [k hop (A) + k hop {B)]/N s (5) 

where khop represents the rate at which each reactant 
migrates from one surface binding site to another, by a 
process of thermal hopping, Ns is the number of bind- 
ing sites on the surface, and f act (AB) is an efficiency 
factor related to the activation energy barrier. This re- 
action rate represents the rate at which two reacting 
species migrate over the surface of a grain and meet 



in the same binding site, and there react. The effi- 
ciency fact(AB) is usually assigned a value equal to the 
probability of reaction upon a single collision between 
the reactants, k(AB), expressed as a Boltzmann factor, 
k(AB) = exp(—EA/T), or a more complex expression 
for quantum tunneling (see section 2.2). However, in re- 
ality, a reaction may have many more opportunities to 
occur, because the reactants are confined in the same 
binding site until one or other reactant hops into an ad- 
jacent site (Tielens & Hagen 1982). Thus, in cases where 
activation energy barriers exist, the competition between 
reaction and migration must be considered. We use an 
expression similar to that given by Chang et al. (2007) 
(different only in that we ignore evaporation, as its rate 
is negligible compared to thermal hopping): 



Jact vk{AB) + k hop (A) + k hop {B) w 

where v is the frequency of collision of the two reactants. 
Here we take v to be equal to the larger of the char- 
acteristic vibrational frequencies of the two reactants, 
as calculated by the method given by Hasegawa et al. 
(1992) and typically on the order of 10 12 Hz. The result 
of equation (6) is that in the limit where thermal hop- 
ping is much faster than the barrier-mediated reaction 
rate, vn{AS), the overall rate of reaction, kAB, depends 
only on the reaction probability, k(AB), and not on the 
hopping rates. Conversely, if the reaction probability is 
sufficiently high, the overall reaction rate is limited only 
by the rate of hopping. 

2.2. Three-phase model 

To the basic chemical model, which consists of a cou- 
pled gas-phase and grain-surface chemistry, we add a 
third phase that represents the ice mantle that forms 
beneath the outer ice surface. This three-phase model 
(gas/surface/mantle) is based on the formulation pre- 
sented by Hasegawa & Herbst (1993). In this formula- 
tion, the surface layer represents the outermost layer of 
ice, which resides either on top of other ice layers, or on 
the surface of the grain itself. The surface and gas-phase 
chemistry arc coupled through the accretion of material 
from the gas phase and the evaporation and non-thermal 
desorption of material from the surface back into the gas 
phase. As the surface coverage of the grain increases, the 
mantle is allowed to grow, incorporating material from 
the surface. The rate of growth of the ice mantle is de- 
termined by the net rate of deposition of material onto 
the grain (which takes into account thermal and non- 
thermal desorption from the grain surface), and by the 
degree of surface coverage of the bare grain. Thus, the 
mantle begins to form even before full surface coverage 
of the grain is achieved, representing deposition of new 
material on top of other surface-bound material. When 
full surface coverage is achieved, any net gain in surface 
material results in a commensurate transfer of material 
to the mantle, with surface material replaced by new 
material from the gas phase. This transfer does not, of 
course, represent a physical motion of material between 
surface and mantle; the transfer is rather an accounting 
issue within the chemical model, as the surface layer is 
continually renewed. 

Following Hasegawa & Herbst (1993), we define the 
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surface concentration of species i as n s (i) = N s (i) ridust, 
where N s (i) is the average, absolute number of parti- 
cles of species i present in the surface layer and ridust is 
the gas-phase concentration of dust particles. Thus, the 
quantity Ng ridust represents one monolayer of ice in this 
parameterization. The calculation of the net loss or gain 
of each species due to chemical and physical processes 
(specifically excluding any loss or gain due to transfer be- 
tween surface and mantle), [dn s (i)/dt] c f iemi is achieved 
using the rates determined in the normal operation of 
the chemical code. The net rate of change in total sur- 
face material is thus [dn s /dt] c hem = ^2i[dn s (i) / 'dt] c hem- 
This quantity is used to calculate both the total rate of 
transfer from surface to mantle, and the rate of transfer 
for individual surface species, respectively, 

[dn s /dt] chem (7) 

[dn s (i)/dt]tran = [dn s /dt] tran [n s (i)/n s ] (8) 

where n s = X)i n s(*), and a acc is the fractional surface 
coverage. The resultant loss/gain rates for the transfer 
of surface/mantle species are appended to the chemical 
rates to give the total loss/gain rates, 

[dn s (i)/dt]tot = [dn s (i)/dt] chem - [dn s (i) / dt] tran (9) 

[dn m (i) / dt]tot = [dn m (i) / dt] chem + [dn s (i)/dt] tran (10) 

where n m (i) represents the concentration of mantle 
species i. 

In the case where [dn s /dt] c hem < 0, and thus there 
is net desorption of material from the ice surface, new 
material is transferred from the mantle to the surface, 
according to the expressions 

[dn s I dt] tran = a des [dn s / dt] chem (11) 



[dn s (i)/dt] tran = [dn s /dt]tran[n m {i)/n m ] (12) 

Contrary to Hasegawa & Herbst (1993), we set ades = 
n m /ns 7^ a acci with ades limited to a maximum value 
of unity. This change allows all surface material to be 
available for desorption, whereas the adoption of otd es = 
o^acc requires the complete removal of the mantle before 
the degree of surface coverage is allowed to fall to less 
than 1. 

In this model we assume the surface material, of a 
thickness of 1 monolayer, to be chemically active, while 
the mantle is taken to be entirely inert. At the grain 
temperatures of around 10 K that are applied to mod- 
els in this study, the mobility of mantle-bound species is 
likely to be negligible. Furthermore, the ability of sur- 
face hydrogen atoms to penetrate the surface and react 
with material in the mantle layers has been shown by 
Fuchs et al. (2009) also to be negligible at such low tem- 
peratures, based on microscopic models of laboratory ice 
experiments. While this treatment is applicable to ices 
detected in quiescent regions, the further chemical pro- 
cessing of layers deep within the ice mantle engendered 
by the star-formation process would require a more com- 
plex chemical treatment of the ice mantle. 



2.3. Activation- Energy Barriers and Quantum 
Tunneling 

It is well known that a number of important grain- 
surface reactions are mediated by activation energy bar- 
riers. The grain-surface formation of CO2, as well as 
that of formaldehyde (H 2 CO) and methanol (CH3OH), 
appears to require reactions with significant activation- 
energy barriers. The O + CO — > C0 2 reaction is fre- 
qcntly invoked as a primary formation route for carbon 
dioxide, due to the availability of the reactants; however, 
it is mediated by a barrier, most recently estimated by 
Roser et al. (2001) to be of height 290 K. 

The activation energies of many pertinent reactions 
have been ascertained with at least moderate precision 
through experiment or calculation. However, a knowl- 
edge of the activation energy alone provides only a very 
weak constraint on the rate of reaction at extremely low 
temperatures. This is because the reaction probabilities 
at around 10 K are so small that quantum tunneling pro- 
cesses can become the dominant means by which a re- 
action proceeds, especially for reactions involving atomic 
hydrogen, and to some degree other atomic species such 
as C and O (see Goumans & Andersson 2010). The rates 
of such quantum processes are generally very poorly con- 
strained. Chemical models such as that of Hasegawa et 
al. (1992) have parametrized the tunneling process by 
adopting an expression for tunneling through a rectan- 
gular potential: 

k = cxp[-2{a/h)(2^E A ) 1/2 ] (13) 

where /x is the reduced mass, and a is the width of the 
barrier, assigned a uniform value of 1 A by Hasegawa et 
al. (1992). Subsequent models have typically taken the 
same approach (e.g. Garrod et al. 2007); however, the 
assumed barrier width is arbitrary. 

Goumans & Andersson (2010) recently calculated the 
temperature-dependent rate of reaction (3) in the gas 
phase, using quantum transition state theory. These cal- 
culations imply a reaction probability per collision at 10 
K of k = 4.8 x 10~ 23 (T. Goumans, private comm.). 
This probability represents a purely quantum tunneling 
process, as the thermal rate is negligible at low tempera- 
tures. Parametrizing this value using equation (13), and 
assuming E A — 290 K (Roser et al. 2001), suggests a 
barrier width of 2.33 A. 

Andersson et al. (2009) calculated gas-phase rates for 
the hydrogen abstraction reaction H + CH 4 — > CH 3 + 
H 2 . That study suggests a reaction probability per col- 
lision at 10 K of k = 3 x 10~ 30 (T. Goumans, private 
comm.). This value is around 15 orders of magnitude 
less than that obtained using equation (13) with a barrier 
width of 1 A and a height of Ea = 5940 K, as assumed 
in previous models. Assuming the same activation en- 
ergy, the Andersson et al. value suggests a square-well 
potential width of 2.17 A. 

Fuchs et al. (2009) conducted laboratory experiments 
on the hydrogenation of CO ice, and modeled this system 
using microscopically exact methods that naturally take 
account of reaction-diffusion competition. They found 
that the H+CO and H+H 2 CO hydrogenation reactions 
proceed by a tunneling process; their combined models 
and data imply a reaction probability of approximately 
k = 10 -14 , to within around 1 order of magnitude. 
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Here we adopt the 10 K tunneling probabilities de- 
rived from the results of Goumans & Andcrsson (2010) 
and Andersson et al. (2009). For the grain-surface hy- 
drogenation of CO and H 2 CO, and in the absence of di- 
rectly measured or calculated rates, we assume the typ- 
ical activation-barrier heights of 2500 K (e.g. Garrod et 
al. 2008), but use a barrier width a = 2 A to calculate the 
square-well tunneling rates according to equation (13). 

2.4. Treatment of the CO +OH^C0 2 +H reaction 

The precise mechanism for reaction (2) involves mul- 
tiple transition states, whose relative energies are some- 
what uncertain (e.g. Goumans et al. 2008, Song et al. 
2006, Frost et al. 1993). However, a simplistic descrip- 
tion of the mechanism is that CO and OH must overcome 
a small activation-energy barrier to form an excited in- 
termediate, HOCO*. The exit channel to form CO2 + H 
from HOCO* is of a similar barrier height to that of the 
reverse entry channel that returns CO + OH. However, 
it is likely that quantum tunneling through this barrier 
makes the formation of CO2 the most probable outcome 
(Frost et al. 1993, Chen & Marcus 2005), if collisional 
de-excitation of HOCO* does not occur first. Goumans 
et al. (2008) suggest that such de-excitation may be pos- 
sible on dust-grain surfaces, through coupling to the ice 
surface, at a rate sufficient to make stable HOCO the 
most likely product. These authors also aver that un- 
der interstellar conditions, reaction of this product with 
atomic hydrogen could lead to a selection of products, in- 
cluding CO2 and formic acid, HCOOH. The large abun- 
dance of interstellar atomic hydrogen, and its high mo- 
bility on grain surfaces, would make such a mechanism 
highly efficient. In such a case, the rate limiting step in 
the formation of CO2 (from CO and OH already present 
in a surface binding site) would thus be the traversing 
of the initial barrier to form HOCO*. However, recent 
experiments by Oba et al. (2010), in which reaction (2) 
was isolated, have indeed shown the CO + OH reaction 
to be capable of forming CO2 on ice surfaces at tempera- 
tures as low as 10 K, while observing no signatures from 
the other suggested product branches. 

Previous chemical models (e.g. Ruffle & Herbst, 2001) 
have assumed an activation energy barrier for reaction 
(2) of 80 K. This appears to be based on the measured 
rate of the gas-phase reaction (Frost et al., 1993), com- 
pared with the collisional rate, giving an efficiency of 
^10~ 4 . Here we adopt the 80 K value, at which the 
thermal reaction probability exceeds the quantum tun- 
neling value obtained through equation (13). In fact, 
the consideration of competition between reaction and 
diffusion described in equation (6) results in a reaction 
efficiency f act that is close to unity in this case, due to the 
relatively high diffusion barrier for CO of 575 K (assum- 
ing Edif = 0.5i?dos)- Thus, in this model, reaction (2) 
is essentially unaffected by the small activation energy 
barrier, and proceeds at very close to the diffusion rate. 
Even adopting the highest values obtained from various 
calculations of the entry channel barrier (e.g. Chen & 
Marcus, 2005), around 400 - 500 K, reaction would prob- 
ably dominate over migration of CO from the binding 
site, especially if quantum tunneling is efficient. 

As well as the production of CO2 by the addition of 
mobile CO to OH, we also incorporate into the model the 
possibility that an oxygen atom may be hydrogenated 



while situated on top of a surface CO molecule. The re- 
sultant chemical energy released would be quite sufficient 
to overcome the small activation barrier to reaction (2), 
and could improve the likelihood of a rapid dissociation 
of HOCO* to C0 2 + H. We formulate this mechanism 
as: 

H + + CO^C0 2 + H (14) 

A similar such mechanism was suggested by Goumans 
& Andersson (2010), based on their calculated gas-phase 
O-CO potential-energy surface, in which the oxygen 
atom may form a loosely bound O...CO complex, allow- 
ing it to remain with the CO molecule for longer, prior 
to the arrival of an H atom. In the case where atomic 
hydrogen is reasonably abundant, such as under typical 
dark-cloud conditions, it is likely that the timescale for 
reaction of H and O would be shorter than that of ther- 
mal diffusion of an O atom. In such a case, the additional 
binding to the CO molecule would be irrelevant, and the 
distribution of O with respect to CO would be essen- 
tially random. Under this assumption, we allow that a 
fraction, #(CO), of all H + O — > OH reactions leads to 
the formation of CO2 + H, where 6>(CO) is the fractional 
surface coverage of CO. This formulation thus represents 
a minimum reaction efficiency for this mechanism, which 
would move (from a value 9(CO)) to something closer to 
unity under hydrogen-poor conditions. 

3. RESULTS 

3.1. Dark cloud ice chemistry 

The method described in section 2 is applied to a static 
single-point dark-cloud model, with density njj = 2 x 10 4 
cm~ 3 and visual extinction Ay = 10. Both gas and 
grains are assigned temperatures of 10 K. Figure 1 shows 
the growth of the ice mantle over time, up to the final 
model time of 10 7 yr. The deposition of the ice is not a 
simple function of time; however, the rate of increase 
over the first 100 monolayers is approximately linear. 
An ice thickness of 1 monolayer is achieved after a few 
thousand years. The rate of mantle deposition falls after 
100 monolayers, as gas-phase material becomes depleted. 
Approximately 300 monolayers have formed after 10 7 yr, 
by which point the (net) deposition rate has fallen close 
to zero. 

Figure 2(a) shows the chemical composition within 
each layer of the ice. A timescale is also shown along 
the top of the plot, which is mapped to the values of the 
lower x-axis according to the relationship shown in fig- 
ure 1. The time corresponding to any particular ice layer 
represents the time in the cloud's chemical evolution at 
which the ice layer was deposited. The chemical com- 
position of the ice beneath the surface layer is preserved 
as the cloud evolves. Thus, the curves to the left of any 
chosen time in the plot represent the instantaneous com- 
position of a section through the ice at that time. 

It may be seen that water is the dominant ice con- 
stituent throughout the model run, and has a relatively 
constant abundance throughout the ice mantle. As found 
by Cuppen & Herbst (2007), the reaction H 2 + OH ->• 
H 2 (E A = 2100 K, Baulch et al., 1984) is the main 
mechanism for the conversion of OH to water, rather 
than reaction with atomic H, because of the large H2 
abundance on the grains (^3 % of surface coverage). 
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Methane (CH4) and ammonia (NH3) are both abundant 
in the first ^50 monolayers of ice, due to the atomic ini- 
tial conditions; C and N are successively hydrogenated 
on the grain surfaces. 

Figure 2(b) shows the gas-phase abundances (with re- 
spect to total hydrogen) of a selection of species, also 
mapped to the ice deposition profile for ease of compar- 
ison. NH3 ice abundance gently falls, through the outer 
layers, as the abundance of CO and its products take up a 
greater proportion of the ice composition. The dramatic 
fall in the abundance of CH 4 ice is caused by the steady 
removal of atomic C from the gas phase, caused primarily 
by its conversion into gas-phase CO. Gas-phase atomic 
oxygen is initially much less affected by its conversion to 
CO, due to its larger abundance. As more CO is formed 
in the gas phase, more is accreted onto the grains, and 
CO reaches ~10 % of the total ice (~20 % with respect 
to water) between monolayers 100 - 250. Gas-phase CO 
reaches its peak abundance at around 3 x 10 5 yr, but it 
only grows to dominate over all species (excluding H 2 ) 
by around 7 x 10 5 yr. 

CO2 ice is formed in appreciable abundances through- 
out the model run, typically 2-4 times less abundant 
than CO. Both CO and CO2 ice abundances follow the 
gas-phase behavior of CO. The only active grain-surface 
CO2 formation route is that of equation (14), in which 
the H + O — > OH reaction occurs on top of a CO 
molecule, leading to the immediate formation of CO2 
+ H. If this mechanism is turned off in the model, CO2 
abundances reach little more than a fraction 10" 3 of CO 
abundance; in that case, reaction (3) is the dominant 
production mechanism. 

The gas-phase abundance of atomic hydrogen stays 
fairly constant at ^10~ 4 n H throughout the model run, 
even as heavier elements deplete from the gas phase. 
However, as the rate of accretion of atomic oxygen from 
the gas-phase falls, the ability of surface oxygen to re- 
move hydrogen atoms on the grains (by the formation of 
OH and H2O) diminishes, and atomic H abundance on 
the grains increases. The fall in water ice production also 
results in a slowing of the rate of ice deposition at late 
times (> 10 6 yr), as seen in figure 1. This reduces the 
rate at which CO (and other surface species) are overlaid 
with new material, leaving them more vulnerable to hy- 
drogenation by accreting H atoms. This combined effect 
of the rise in grain-surface atomic hydrogen abundance 
and the slower rate of ice deposition leads to a much 
greater efficiency in the hydrogenation of CO and H 2 CO 
to methanol (CH3OH). The effect is further enhanced as 
CO takes over as the dominant accreting species other 
than hydrogen, at ~1 Myr. Thus, the formation of large 
quantities of methanol in the outer ice mantle layers is 
a symptom of the depletion of gas-phase oxygen and its 
conversion in the gas-phase into CO, under conditions 
where atomic H is abundant. Charnley & Rodgers (2009) 
found a similar effect in their model, albeit apparently in- 
fluenced by the action of hydrogen abstraction reactions 
on surface methanol and its precursors. 

Figure 2(b) also shows the abundance of gas-phase 
methanol, function of the number of ice layers in 
the mantle. The CH3OH fractional abundance matches 
the canonical dark-cloud value of ~10~ 9 , from around 
10 5 yr until near to the end of the run. This methanol is 
formed on the dust-grain surfaces, and is ejected into the 



gas phase with an efficiency of < 1% by the chemical en- 
ergy released upon its formation, under the formulation 
given by Garrod et al. (2007). At the 10 magnitudes of 
visual extinction used in this model, direct photodesorp- 
tion at the yields determined by Obcrg et al. (2009a,b) is 
found to be negligible compared to reactive desorption. 
While the conversion of CO to CH 3 OH on the grains be- 
comes more efficient at later times, the overall rate of 
conversion remains fairly steady, due to the falling rate 
of deposition of CO onto the grains from the gas phase. 
This leaves the gas-phase abundance of methanol stable 
over a long period. 

Molecular oxygen is formed in the gas-phase, mainly 
by the reaction O + OH — > O2 + H; molecular oxygen 
abundance is shown in figure 2(b). This 2 can accrete 
onto the grain surface and contribute to the formation 
of surface OH and water. This occurs primarily by the 
sequential hydrogenation of 2 by atomic H, to produce 
hydrogen peroxide, H 2 2 , which is followed by the hy- 
drogen abstraction reaction H + H2O2 — > H2O + OH. 
The reactions of H with O2 and H2O2 are assigned ac- 
tivation energy barriers of 1200 K and 1400 K respec- 
tively (Tielens & Hagen 1982), following the networks of 
Ruffle & Herbst (2001) and earlier models by those au- 
thors. However, the treatment of competition between 
reaction and diffusion used in the present model (see sec- 
tion 2.1), combined with quantum tunneling through the 
activation barriers (section 2.3), allows both reactions to 
proceed at close to the hydrogen diffusion rate. This 
high efficiency is in line with experimental findings by 
Miyauchi et al. (2008), Ioppolo et al. (2008), and Oba 
et al. (2009) showing that O2 may be efficiently con- 
verted to water at low temperatures. In this model, the 
O2 to H2O conversion is sufficiently rapid that there is 
no significant build-up of O2 (or indeed H2O2) in the 
ice mantles. The contribution of 2 -related chemistry 
to the water-ice formation rate is around 1 % until close 
to the end of the model run, when gas-phase O2 briefly 
dominates over atomic O. At its peak, the H + H2O2 — > 
H 2 + OH reaction contributes around 30 % of water 
production, and ~ 50 % of OH production. The major 
O2 contribution to water formation is thus confined to 
the outermost ice layers, and has no discernible effect on 
the water profile through the ice, as shown by figure 2(a). 

3.1.1. Variations in physical/ chemical parameters 

Here we test the effect on the basic dark-cloud chemical 
model of variations in the dust temperature, the cloud 
density, and the generic diffusion-barrier/binding-encrgy 
ratio. 

Figure 3 shows the ice mantle composition for model 
runs with dust temperatures of 8 K and 12 K. The gas 
temperature is held constant at 10 K, in order to isolate 
the surface/ice-specific effects. In fact, the general be- 
havior of the chemistry in the 8 - 12 K range is quite ro- 
bust. The main difference between models, as judged by 
the observationally-dctectable, saturated species shown 
in the figure, is that the efficiency of methanol production 
is generally lower at higher temperatures. This is due to 
the greater mobility of hydrogen, which increases its reac- 
tion rate with surface radicals and thus lowers its surface 
abundance, while leaving unaffected the total reaction 
rate for barrier-mediated reactions such as H + CO — > 
HCO and H + H 2 CO -> CH 3 0/CH 2 OH. Because both 
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the formation and destruction of formaldehyde (H2CO) 
are governed by these reactions, the abundance of H2CO 
relative to CO is less affected by temperature changes, 
and broadly follows the behavior of CO. 

The constancy of general behavior over the 8 - 12 K 
range shown in figure 3 is largely due to the system being 
in the accretion limit. This means that surface reactions 
(without activation energy barriers) are rapid compared 
to the rate of accretion, and thus production occurs at 
these limiting accretion rates rather than according to 
the temperature-dependent surface-diffusion rates. How- 
ever, the averaged abundances of intermediate species 
such as OH or CH 3 are strongly affected by the tem- 
perature variation; their small surface abundances (i.e. 
typically less than one particle per grain at any instant) 
reflect their lifetimes on the grain surface before reac- 
tion occurs. Reactions with activation energy barriers 
(which are typically much slower than accretion rates) 
are also largely unaffected by temperature changes, due 
to the dependence of their rates primarily on the acti- 
vation barrier tunneling rate, as per equation (6). The 
variation with temperature of the overall number of ice 
layers formed is caused by the larger proportion of H 2 
being incorporated into the ice mantles at lower temper- 
atures, due to the lower rate of evaporation. 

Figure 4(a) shows the results of the model using a 
generic assignment of £?<jif = 0.3-Edos ( as adopted by 
Hasegawa et al. 1992, based on estimates given by Tie- 
lens & Allamandola 1987) with a dust temperature of 
8 K. Again, the change in the general chemical behav- 
ior is small; the lowering of diffusion barrier heights is 
analogous to an increase in dust temperature, in the 
regime modeled here. The main difference in behavior 
compared with the E dii = 0.5i?dcs case, is the initially 
higher methanol ice abundance, which is unrelated to 
CO surface abundance. This methanol formation is en- 
abled by the increased mobility of oxygen atoms, which 
react at a sufficiently fast rate with methyl (CH3) radi- 
cals to contribute significantly to CH 3 formation, and 
thus to methanol production. A similar increase in sur- 
face mobilities is found in the basic model with a dust 
temperature of 12 K, as seen in figure 3(b); however, 
in that case, the abundance of CH3 is lower due to the 
much lower molecular hydrogen abundance on the grains, 
whose presence otherwise provides a mechanism for the 
hydrogenation of C to CH 2 and of CH 2 to CH 3 . This 
alternative methanol-production route falls away as the 
rate of accretion of atomic carbon onto the grains dimin- 
ishes. 

Figure 4(b) shows the same model (E^f = 0.3i?dcs), ex- 
cept with an increased dust temperature of 12 K. In this 
case, the chemical behavior on the grains is far different; 
CO2 dominates CO throughout the mantle, and even be- 
comes the dominant ice-mantle constituent through most 
of the ice layers. In this case, CO has become sufficiently 
mobile for reaction (2) to compete with H 2 for reaction 
with OH; the elevated temperature also decreases H 2 cov- 
erage, which stands at less than 0.01 % at 12 K. 

Clearly, under certain temperature conditions, or with 
somewhat higher mobilities of surface species, C0 2 can 
be seen to dominate the ice composition. However, at 
no point in the model shown in figure 4(b) does CO take 
over from C0 2 ; the use of a single-point static chemi- 
cal model appears to be insufficient to obtain a strongly 



time-dependent relationship between CO and C0 2 abun- 
dance in the ice. 

Figure 5 shows ice-mantle and gas-phase abundances 
using the basic £dif = 0.5i?des model with a density re- 
duced by a factor of ten, to rin = 2 x 10 3 cm~ 3 . Here 
both CO and C0 2 comprise only a few percent of the 
total ice mantle, while hydrogenation of CO to methanol 
is highly efficient. This is caused by the ten times higher 
atomic-hydrogen abundance that results from the ten 
times slower conversion of H to H 2 . Here, accretion of 
atomic H onto the grains is much faster than that of 
atomic O, thus there is little competition for hydrogen 
between O and CO, and surface CO is rapidly converted 
to methanol at all times. In this low density case, the 
depletion of heavy elements does not reach completion 
by 10 7 yr. The somewhat lower rate of mantle formation 
also allows more time for CO hydrogenation before the 
CO molecules are covered by new ice layers. 

Figure 6 shows ice-mantle and gas-phase abundances, 
assuming a density n H = 2 x 10 5 cm~ 3 . Here, the gen- 
eral behavior of the chemistry is more similar to the 
«h = 2 x 10 4 cm~ 3 case, with the exception of methanol. 
Atomic hydrogen abundance in the gas phase is typically 
less than that of heavier elements or CO. Thus, oxygen 
reacts with most of the surface H, leaving the majority 
of CO intact. Only at very late times - in the outermost 
ice layers - does methanol become abundant, as heav- 
ier species become depleted and atomic hydrogen again 
dominates the accretion from the gas phase. Gas-phase 
methanol, which originates on the grains, never reaches 
the typical observed value of ^10~ 9 nn- Total depletion 
of most gas-phase species occurs at around 10 6 yr. 

3.1.2. Discussion of dark cloud model results 

The static dark-cloud models show that a significant 
fraction of CO may be converted to C0 2 on the dust 
grains at low temperatures. This behavior appears quite 
robust over a moderate range of dust temperatures and 
diffusion barrier heights. The mechanism at work is re- 
action (2); however, the means by which it occurs is not 
through the diffusion of the reactants CO and OH, but 
rather the formation of OH in situ in a binding site that 
already contains a CO molecule. Thus, the only surface 
diffusion necessary for the formation of C0 2 in this case 
is atomic hydrogen; atomic O may simply remain where 
it lands following accretion. The proportion of CO that 
is converted to C0 2 is then simply dependent on the 
fraction of the surface composed of CO. 

The suggestion by Goumans & Andersson (2010), that 
O and CO may also form a loose bond, could increase 
the conversion rate, where the atomic hydrogen accre- 
tion rate is low enough to allow atomic oxygen time to 
diffuse significantly around the grain surface. We also 
tested this mechanism using the basic 10 K dark cloud 
model, adding a new species, O...CO, to the model to 
represent this loosely bound association of O and CO, 
and allowing reaction with a selection of prevalent re- 
action partners including atomic H. However, we found 
that this approach yielded C0 2 production efficiencies 
lower than those obtained with the method outlined in 
section 2.4. This was due to the lack of an Eley-Rideal 
mechanism for the formation of O...CO, whereby oxygen 
could accrete directly onto a CO molecule - the forma- 
tion of O...CO in the model was thus solely dependent 



8 



Garrod & Pauly 



on the diffusion of atomic O. In the case where such an 
Eley-Rideal mechanism is the dominant O...CO produc- 
tion route, the result should reduce to that of the sim- 
pler approach in any case. The incorporation of O...CO 
into the reaction network in a sufficiently comprehensive 
fashion, including an appropriate treatment under the 
modified-rate method of Garrod (2008), is left for future 
exploration. 

At somewhat higher temperatures and/or with lower 
diffusion barriers, the efficiency of CO to CO2 conver- 
sion undergoes an enormous increase, as CO becomes 
sufficiently mobile on the grain surfaces to compete for 
reaction with OH - reaction (2) becomes active through 
direct surface diffusion. It is therefore plausible that 
higher temperatures during the early formation of the 
ice mantles could produce the highly efficient C0 2 for- 
mation observed at low visual extinctions. 

The ice composition produced by the models is in fair 
agreement with observational values. Table 3 shows the 
ice observations, as collated by Gibb et al. (2000), to- 
ward the field star Elias 16, which is commonly taken as 
a representative dark cloud source. Alongside, we give 
the calculated composition obtained with the basic 10 
K model, and summed over all the ice layers present at 
a given time, for a range of model times. Of the times 
shown, 10 6 yr might be chosen as the optimum; however, 
the abundances of CH4 and NH 3 are excessively high at 
all times in the model. Judging by the results shown in 
figure 4, this could perhaps be remedied somewhat by the 
choice of a lower dust temperature and diffusion barrier 
height. Such changes would have no significant effect on 
the agreement of CO/CO2 abundances, as seen in figure 
4(a). The higher temperature (12 K) case seen in figure 
4(b) demonstrates an improved CH 4 agreement, but NH 3 
agreement is generally worse, while the CO/CO2 behav- 
ior shows strong CO2 dominance as mentioned above, 
which does not fit with observational results. However, 
since both CH4 and NH3 are particularly abundant early 
in the model cloud's evolution, it is also likely that a bet- 
ter dynamical treatment of the cloud's evolution at early 
times would improve the fidelity of the model to the ob- 
served values. Both CH4 and NH3 abundances within the 
ice are lower (i.e. better), under the 12 K regime, at early 
times (or for the deepest ice layers), while lower tem- 
peratures show better agreement at intermediate times 
(middle ice layers) . This would suggest that the collapse 
and cooling of the cloud may play a significant role in 
controlling NH 3 and CH4 abundances in the ice. 

Methanol and formaldehyde abundances are also a lit- 
tle high at later times, comparing the basic 10 K model 
with the observational values, but variation from the as- 
sumed cloud density could change this. 

3.2. Depth- dependent models 

Here we attempt to model the visual extinction- 
dependent chemistry of water, CO and CO2 ices, in order 
to address the observed extinction-threshold behavior of 
these species. We assume a constant density through- 
out a dark cloud, and apply the model at a selection 
of depth points from the edge of the cloud toward the 
center, corresponding to fixed visual extinctions. Cru- 
cially, we adopt a visual extinction-dependent dust tem- 
perature, assuming exposure to the standard interstellar 
radiation field (ISRF). To calculate this temperature in 



a simple manner, we assume that the rate of cooling of 
the dust is equal to the rate of radiative heating by the 
ISRF, thus we solve: 



QMT dust )dv= / Q v J v D v (Ay)dv (15) 
Jo 

where Q v is the frequency-dependent efficiency of emis- 
sion and absorption of the grain, -B„(Td us t) is the Planck 
function for dust-grain temperature Tdust, Jv is the ra- 
diation field intensity, and D v {Ay) is the attenuation of 
the field at given frequency by the dust column ahead 
of the position in question. To simplify this calcula- 
tion, we evaluate the left-hand side using the approxi- 
mation of Kriigel (2008), L.H.S. ~ 1.47 x 10~ 6 aT$ uat , 
expressed in cgs units. To evaluate the right-hand side, 
we use the radiation field of Zucconi et al. (2001), which 
consists of multiple, modified black-bodies ranging from 
the optical-NIR down to cosmic microwave background 
radiation. To evaluate the emission efficiency we use 
Q v ~ 10~ 23 av 2 (Kriigel 2008), with v in Hz and a in cm. 
We follow the approach given by Cuppen et al. (2006) 
to evaluate the attenuation of the radiation field, 



D u = exp 



-0.8 



Ay 



(16) 



using the tabulated values given by Mathis et al. (1990) 
to evaluate A v jAy. Our treatment implicitly assumes 
plane-parallel geometry. Figure 7 shows the resulting 
dust temperature profile, using the canonical 0.1 [im 
grain size; Td us t = 8.41 K at Ay = 10. 

At low visual extinctions, the photodissociation rates 
of H2 and CO become important to the chemistry. We 
use the tabulated values of Lee et al. (1998), which de- 
pend on visual extinction and on the column densities of 
H2 and CO. Beginning at the edge of the cloud, the code 
works inward, summing the column densities through 
previously calculated points (at a given time step), to 
give the current column density values. The final ratio 
of H:H 2 for each depth point is used as the starting ratio 
for the next point inward. The model was run over a 
visual extinction range of - 10 magnitudes, from cloud 
edge to center. A linear relationship between total hy- 
drogen column density and visual extinction is assumed, 
with iV(H+2H 2 ) = 1.6 x 10 21 Ay. We adopt a fixed gas 
temperature of 10 K. 

Figure 8 shows the depth-dependent column-density 
behavior of H2O, CO2 and CO ices for three constant 
gas densities: iir = 2x 10 3 (low), 6 x 10 3 (intermediate), 
2 x 10 4 cm~ 3 (high), assuming a generic grain-surface 
diffusion-barrier height -Eyif = 0.3£y C s- Curves for each 
molecule are shown for three different times in the mod- 
els: 0.5, 1 and 2 Myr; higher curves correspond to later 
times. The visual extinction values shown in the figure 
are twice those used in the chemistry model, so as to rep- 
resent edge-to-edge values. The column density at any 
given visual extinction represents the sum of the calcu- 
lated abundances, from the edge of the cloud inward to 
that point, multiplied by 2 to give an edge-to-edge value. 
Thus, the plotted column densities and visual extinctions 
are directly comparable to the observed quantities. The 
dotted line represents an order-of-magnitude estimate of 
the detection threshold for all species. More precise es- 
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timates may be found in the literature (as cited in the 
Introduction) . 

It is immediately seen that the general, observed 
visual-extinction threshold behavior is reproduced by the 
models; water ice is formed at the lowest visual extinc- 
tion of each of the three molecules, with CO2 reaching the 
threshold (10 16 cm -2 ) at marginally higher Ay, and CO 
following at significantly greater visual extinctions. The 
threshold for H 2 in the low- and intermediate-density 
runs is a remarkably good match to the observational 
value, with values of 4 and 3 magnitudes, respectively. 
The water threshold is fairly insensitive to the age of the 
cloud, over the 0.5-2 Myr range shown. Water is the 
dominant ice constituent, and its threshold visual extinc- 
tion represents the onset of the deposition of more than 
1 layer of ice on the grains. At lower extinctions, the bal- 
ance between desorption and accretion holds the surface 
coverage of the grain to less than unity. Conversion of 
CO to CO2 is efficient at low visual extinctions, due to 
the higher dust temperatures, which allow CO to diffuse 
sufficiently rapidly to react with surface OH via reaction 
(2). For the high density run, the threshold visual ex- 
tinction for H 2 is significantly lower, at around 1.5. We 
therefore suggest that the water (and CO2) threshold vi- 
sual extinctions are representative of material with a gas 
density of a few times 10 3 cm~ 3 . 

Figure 9 shows column densities calculated for mod- 
els with various Edu'.Edca ratios, adopting a gas den- 
sity of 6 x 10 3 cm~ 3 . Panel (a) shows results for the 
Edif — 0.35 E dcs case. Comparing with figure 8(b), we 
see that the slightly higher diffusion barrier moves the 
CO extinction threshold to lower values, while those of 
H2O and CO2 are essentially unchanged. A further in- 
crease to E d [f — 0.4E dos , shown in panel (b), results in 
a CO2 threshold in good agreement with the observed 
value of 6.7 ± 1.6 (Whittet et al. 2007), with the precise 
value dependent on the cloud age. Panels (c) and (d) 
show the results for Edif:Ed es ratios of 0.45 and 0.5, re- 
spectively. Here, the behavior of CO and CO2 changes; 
the lower mobility of CO in these cases renders inefficient 
the formation of C0 2 by direct diffusion and reaction 
with OH, leaving the main CO2 formation mechanism as 
that in which OH forms in a binding site where CO is 
already present (equation 14). In these cases, the thresh- 
old of CO lies at a lower value than that of C02- Thus, 
only values E^if < 0.4 Ed os are consistent with the obser- 
vational data. However, higher visual extinctions may in 
reality have higher densities than do the edges, meaning 
earlier onset of CO formation and a lower threshold ex- 
tinction than calculated here (i.e. closer to the observed 
value). 

For easier comparison with observational data, linear 
plots of the visual extinction-dependent column densi- 
ties for the cases Edif = 0.3Ed os and E^a — 0.4Ed os 
are shown in figure 10. While the curves for H 2 and 
CO ices are fairly linear, it may be seen that the CO2 
curves become less steep at higher visual extinctions, 
contrary to observational data (e.g. Whittet et al., 
2007). This occurs at around 10 magnitudes of extinc- 
tion (edgc-to-edge) , where appreciable CO abundances 
are found. This effect is explainable by the lack of dy- 
namical considerations in the models. The visual ex- 
tinction for a parcel of gas is fixed throughout its evolu- 
tion, so that high visual-extinction positions never expe- 



rience the low-extinction/higher-temperature conditions 
required to form CO2 efficiently. The inclusion of dy- 
namics in these models would likely improve the behav- 
ior at high visual extinctions, by allowing a gradual in- 
crease from low to high extinction and from high to low 
dust temperatures. The C0 2 formed at lower extinctions 
would be incorporated and preserved in the deepest ice 
mantles, allowing it to contribute to the total CO2 col- 
umn through the cloud later on. 

While the intercepts of the linear plots for both CO2 
and H 2 are very similar for all cloud ages in the 0.5 - 
2 Myr range that is plotted in figure 10, the gradients of 
these curves vary with time. Whittet et al. (2007) obtain 
a slope for CO2 of 0.252 x 10 17 cm~ 2 mag -1 , using data 
for various lines of sight in the Taurus dark cloud. Using 
the E d [{ = 0.4E des model shown in figure 10(b), which 
has approximately the correct intercept for CO, the CO2 
gradient of the model over the 3 < Ay < 8 range best 
matches the observed value at a cloud age of 1.1 Myr. 
A similar comparison using the H 2 slope of 1.5 x 10 17 
cm~ 2 mag -1 for the Taurus dark cloud given by Chiar 
et al. (1995) indicates a cloud age of 1.6 Myr. Chiar et 
al. also give H 2 column density data for the Serpens 
cloud. Using the same simple comparison to our model, 
we estimate a cloud age of 1 Myr, somewhat younger 
than either estimate for Taurus. Of course, without con- 
structing a specific physical model for either cloud, such 
estimates are rather imprecise. 

3.3. Collapsing cores 

In order to understand the chemical behavior of the 
ice under collapse conditions, in which visual extinction 
grows and thus dust temperature falls, we employ the 
simple free-fall collapse model given by Brown et al. 
(1988) and based on the construction of Spitzer (1978). 
We begin collapse at a gas density iih,o = 3 x 10 3 cm~ 3 , 
and with a visual extinction of Ay$ = 2. A peak den- 
sity is chosen at which, once reached, the model remains 
for the duration of the run. The visual extinction evo- 
lution obeys the relationship Ay = ^4v,o(*ih/"-h,o) 2 ^ 3 - 
Our adoption of a free-fall collapse model with an arbi- 
trary maximum density is neither expected nor intended 
to represent the complexity of cloud-collapse dynamics, 
but provides a convenient treatment to understand the 
general behavior of the cloud chemistry under such con- 
ditions. 

To ensure reasonable accuracy in the calculated dust 
temperatures at high visual extinctions, we use the dust- 
temperature profile calculated by Zucconi et al. (2001) 
(their equations 22-24 and 26) for extinctions Ay > 10, 
and the profile derived above in section 3.2 for extinctions 
Ay < 10. In order to allow the code easily to sample from 
our profile at every time step in the calculation, we fit to 
the Ay < 10 profile a third-order polynomial: 

T dust = 18.67- 1.637 A v + 0.07518 - 0.001492 

(17) 

To this profile we add a corrective of 0.316 K to en- 
sure continuity at Ay = 10. Zucconi et al. suggest their 
profile is valid in the Ay = 10 - 400 range (but they do 
not define temperatures between 2 and 10 magnitudes of 
extinction, hence the use of our own approximate pro- 
file). We again use the H 2 and CO photodissociation 
treatments given by Lee et al. (1998). 
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Figure 11(a) shows ice mantle composition for collapse 
to a density of 4x 10 4 cm -3 , which wc identify with a dark 
cloud. This density is chosen to produce a visual extinc- 
tion of approximately 10 magnitudes, which results in a 
dust temperature of ~8.7 K. A diffusion barrier height 
of 0.3£?dos is chosen. One immediately sees that C0 2 is 
formed in great abundance at early times, building up 
around 50 layers in a mixture mainly composed of H2O 
and CO2; CO abundance is minimal in this regime. CH4 
also reaches its peak fractional composition in these lay- 
ers, as does ammonia. CO2 reaches its peak at around 
7 x 10 5 yr, at which the dust temperature is 12 K. After 
this, there is a rapid switch-over in behavior, and CO 
dominates over CO2. We identify the CO2 formed in 
layers 0-50 with the observed polar CO2 component. 
At temperatures of less than 12 K, CO diffusion becomes 
unable to compete with hydrogenation of OH, and so the 
alternative mechanism of OH formation in the presence 
of CO becomes the most important for CO2 production. 

Figure 11(b) shows the gas-phase abundances for the 
same model, also mapped to the number of ice layers. 
Gas-phase atomic carbon is seen to fall off at early times, 
as it is converted to CO at relatively low gas densities. 
Initially high gas-phase carbon abundances produce the 
early CH4 ice peak. This declines with increasing rapid- 
ity as the gas density rises, accelerating the conversion of 
gas-phase C to CO. This coincides with the sharp fall in 
CO2, also caused by rising density and visual extinction, 
and thus dust temperature. These two species should 
be well correlated in ice observations. This behavior 
differs from the static dark-cloud models, which forced 
gas-phase CO to form concurrently with the ice man- 
tles, producing high methane content in the ice. Under 
collapse conditions, much atomic carbon is converted to 
CO before significant ice build-up occurs. Even longer 
periods at low density/ Ay prior to collapse would allow 
more carbon to be locked into CO, resulting in yet lower 
methane ice abundances. 

It is also seen that methanol, formed on the grains, 
maintains an appropriate gas-phase abundance, in the 
region where CO dominates C0 2 in the ice mantles. 
Therefore, in regions of visual extinction less than the 
CO threshold, no appreciable gas-phase methanol is to 
be expected. 

Figure 11(c) shows the cumulative fractional composi- 
tion of the ice summed through the entire ice mantle (up 
to a thickness of the value shown on the abscissa), with 
respect to water. These values represent only the ice 
composition within a single parcel of cloud gas, and do 
not take into account any integration through the outer 
parts of the cloud. For an ice thickness of around 50 
layers, cumulative CO2 fraction is at its peak, compris- 
ing approximately 80 % of the water abundance. This 
value falls for greater ice thicknesses, eventually falling 
below the total CO fraction. The last column of Table 
3 shows the "best-match" results for this model as com- 
pared to the dust-grain ice composition for Elias 16 pre- 
sented by Gibb et al. (2000). We roughly determine the 
"best-match" as the time at which both CO and CO2 
comprise 30 % of the ice mantle, which occurs for ice 
thickness of ^220 layers. The new values are a much 
better match, even for CH4, which was problematic for 
the static dark-cloud models. The abundance of NH 3 is 
nevertheless rather high for dark-cloud observational val- 



ues, although it does not exceed levels observed toward 
star-forming regions (20 - 30 %, Gibb et al., 2000), after 
its initial peak. The introduction of collapse from a lower 
visual extinction also allows a much better reproduction 
of observed CO2 values. The values determined by the 
model are of course somewhat dependent on the choice 
and complexity of the dynamical model. 

Figure 12 shows the results for collapse to a density of 
10 7 cm -3 , which is achieved after around 9.5x 10 5 yr, and 
by which time extreme depletion of gas-phase species has 
occured. The general composition and behavior of the ice 
is not greatly different from the lower density run; how- 
ever, methanol formation falls dramatically after its peak 
at around 50 layers. This is caused by the much lower 
dust temperature (~ 5 K), which allows H2 to maintain 
a dominant level of ice-surface coverage. The high depo- 
sition rate of H2 onto the ice reduces the period of time 
available for atomic hydrogen to hydrogenate CO and 
H2CO to methanol before a new layer of ice is deposited 
on top. Thus, the formation of significant quantities of 
methanol in very dense cores, with such low tempera- 
tures, is unlikely. However, it is not clear whether dust 
temperatures this low could be achieved before star for- 
mation were to warm the dust grains, and the free-fall 
collapse model produces only a simplistic treatment of 
physical conditions. 

4. DISCUSSION 

These models show that CO2 may be formed on dust- 
grain surfaces by the CO + OH — > CO2 + H reaction 
alone. The efficiency of this barrier-mediated reaction is 
determined by the competition between reaction and the 
loss of the reactants by their diffusion out of the surface 
binding site. Using our standard surface-binding char- 
acteristics, we find the efficiency to be close to unity. 
Under such a regime, the specific mechanism and rate 
at which CO and OH meet on the surface are the de- 
termining factors for CO to CO2 conversion. At suffi- 
ciently high temperatures (~12 K), CO is mobile enough 
to reach and react with OH radicals before either H or 
H 2 , producing a high CO-to-C0 2 conversion efficiency 
and allowing CO2 to dominate CO in the ice mantles. 
At lower temperatures, CO2 is still formed in significant 
quantities; however, in this case the OH is formed atop a 
CO molecule by the addition of O and H, and CO abun- 
dance always dominates that of CO2. Typical gas-grain 
codes whose reaction rates are determined by diffusive 
processes do not take such a mechanism into account. 
Importantly, the only diffusion required for this process 
is that of atomic hydrogen; oxygen may be hydrogenated 
wherever it lands on the grain surface. While the release 
of chemical energy associated with the formation of OH 
is likely to render the CO + OH reaction more rapid 
(perhaps increasing the probability of the dissociation of 
the excited intermediate HOCO* into C0 2 and H), in 
this model it is simply the act of bringing CO and OH 
together that allows the mechanism to be so efficient. 

In contrast to reaction (2), reaction (3) is rendered 
ineffective by the calculated quantum tunneling rates 
of Goumans & Andersson (2010), while reaction (1) is 
never competitive, due to the relatively low mobilities 
and abundances of each reactant (O and HCO). 

It is of note that time-dependent changes in the gas- 
phase chemistry of static dark-cloud models do not ap- 
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pear able to induce a switch-over from CO2 to CO dom- 
inance in the ice mantles; it would seem that a change 
in physical conditions is required for the CO:C02 ratio 
in the ice mantles to be drastically altered. While the 
gas density does not appear to have a direct effect on 
the CO:C0 2 ratio, adjustments to the dust temperature 
are capable of inducing a switch-over, by changing the 
dominant mechanism by which CO and OH meet on the 
grains, as mentioned above. However, at temperatures 
in the 8 - 12 K range we find that the general behavior 
of the ice chemistry is remarkably robust; previous mod- 
els (e.g. Ruffle & Herbst 2001) have often shown strong 
variation with temperature, due to the use of the sim- 
ple deterministic "rate-equation" approach to solve the 
surface chemical network. 

Our simple treatment of the dependence of dust tem- 
perature on visual extinction allows the observational 
threshold extinctions of H2O and CO2 to be well re- 
produced, assuming a gas density of a few xlO 3 cm~ 3 . 
Higher densities tend to lower the threshold extinction 
for these species. The threshold extinction for CO can 
also be reproduced by the models, but the precise value 
is dependent on the ratio of diffusion rates to binding en- 
ergies on the grain/ice surfaces, and is also more sensitive 
to the assumed age of the cloud than those for the other 
two molecules. The diffusion barriers determine the tem- 
perature at which surface CO becomes sufficiently mobile 
to induce a switch over in CO:C02 dominance. However, 
using the adopted dust temperature profile, the barriers 
must be less than 0AE des to be able to produce the 
switch-over at all. Somewhat lower values than this can- 
not be ruled out, as the use of a more realistic physical 
model - either through collapse or through a non-uniform 
density profile - may also cause the modeled CO:C0 2 
boundary to shift. 

The observational extinction threshold for H2O corre- 
sponds to the point where the photodesorption of water 
(and other surface species) back into the gas phase is not 
sufficient to hold the ice to less than one layer of surface 
coverage. As the rate of ice formation increases further, 
the total photodesorption rate is held steady (in absolute 
terms), since material can only desorb from the surface 
and not from deeper layers. This effect will remain es- 
sentially the same even if photodesorption occurs from 
the first few surface layers. 

C0 2 formation occurs as soon as gas-phase CO is 
present. In our low-density collapse models, a few lay- 
ers of water ice are formed before gas-phase CO becomes 
abundant. However, a collapse model that began at a 
visual extinction lower than 2 magnitudes would show 
CO formation before any significant ice mantles were 
formed. The slightly higher observed extinction thresh- 
old for C0 2 than for H2O probably derives from the 
slightly lower abundance of CO2 in the ices. 

In the case of CO, the extinction threshold is tied to 
a temperature threshold of around 12 K, at which point 
CO mobility is insufficient for reaction (2) to occur via 
direct diffusion. As a collapsing cloud or core reaches 
its threshold visual extinction, the dust becomes cool 
enough that the less efficient mechanism takes over. 

The efficient formation of CO2 by the addition of CO 
and OH via direct surface diffusion is dependent on the 
reactants having sufficient time in the binding site to re- 
act before one of them diffuses away. If the barrier to 



the first stage of this reaction, corresponding to CO + 
OH — > HOCO*, is low enough enough then formation of 
CO2 may occur. There is some debate as to the precise 
value of this barrier, however it is typically estimated to 
be low enough that our assumption of high efficiency for 
reaction (2) is very plausible. As discussed by Goumans 
et al. (2008), the rapidity with which excited HOCO* 
loses energy to the grain surface may be important to 
the mechanism by which CO2 is ultimately formed. But, 
the abundance of gas-phase H combined with its rapid 
surface diffusion make quite probable the alternative pro- 
cess of H-addition to de-excited HOCO that those au- 
thors suggest. A more detailed chemical model would 
be valuable in confirming this expectation. Goumans et 
al. (2008) also suggested that HCOOH would be another 
significant product branch for HOCO hydrogenation. No 
alternative product branches are included in our network. 
However, Oba et al. (2010) detected no such signatures, 
perhaps indicating that the loss of energy to the surface 
is not sufficiently fast to affect the direct dissociation of 
HOCO to H + CO2. All of these considerations are im- 
portant in the case of direct diffusion of CO to OH, but 
the mechanism by which O and H combine on top of a 
CO molecule would involve a highly excited OH radical 
which could overcome both any initial barrier to HOCO* 
formation and the internal barrier to the ejection of an 
H atom to produce CO2. 

The formation of CO2 at low visual extinctions, as 
demonstrated by the models, requires a certain max- 
imum diffusion-barrier height of E di{ ~ 0AE dcs for 
CO, although here we change the generic Edif.Edcs ra- 
tio to avoid excessive parameter space. It is possible 
that higher values could produce the same effect, under 
higher dust-temperature conditions, although our simple 
extinction-dependent temperature profile agrees reason- 
ably well with the much more precise approach of Zucconi 
et al. (2001) whose data we use to make our approxima- 
tion. However, those authors do not give well-defined 
temperatures in the Ay=2 - 10 range. A departure from 
the profile we assume could thus affect the threshold ex- 
tinction for CO produced by the models. 

CH4 is generally less abundant than NH3 in the ices 
due to the gas-phase conversion of atomic carbon to CO, 
where atomic nitrogen does not have a similarly rapid 
removal process. Water production is not affected in the 
same way as CH4, due to the larger gas-phase abundance 
of atomic oxygen. In the case of a collapsing core, CH4 
and C0 2 are correlated in the ices (a result also suggested 
by the observational data of Oberg et al., 2008), due to 
the formation of CH4 at early times, and the formation of 
C0 2 at the high temperatures and low visual extinctions 
of the pre-collapse stage. 

Observations indicate that a large proportion of in- 
terstellar CO2 ice is embedded in a polar mixture (e.g. 
Gibb et al. 2000, Bergin et al. 2005), i.e. a water-rich 
ice, and this component of the ice mantles may be eas- 
ily explained by the efficient conversion of CO to CO2 
at low visual extinctions and elevated temperatures. In 
dark clouds, this component would thus be formed at dy- 
namically early times, before the gas density and visual 
extinction had reached more typical dark-cloud values. 
Furthermore, observations toward the Taurus field star 
Elias 16 by Bergin et al. (2005) found that, while the ma- 
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jority of detected CO2 ice is in a polar mixture, around 15 
% resides in a non-polar CO:C02 ice dominated by CO; 
indeed, their best match to this ice composition was a 
100:26 mixture. The static dark-cloud models we present 
here typically produce a CO:C02 ratio of approximately 
4:1 to 2:1, dependent on the physical parameters of the 
model - this ratio is very close indeed to the (non-polar) 
observational result, and is fairly robust through many 
layers of ice. However, at all times in those model runs, 
water is the most abundant ice species. Indeed, there is 
no obvious means by which one can engineer the models 
to produce both a low water content and an appropriate 
CO:C02 ratio that endures for a sufHcent time. Under 
the dominant CO2 formation mechanism of the dark- 
cloud models, in which OH is formed in the presence of 
CO, a significant amount of water is always formed, as 
the price of the inefficiency of CO to CO2 conversion 
that produces the CO:C02 ratio of ~2 - 4. It is possible 
that mobile atomic oxygen could become more strongly 
bound to CO, such that they remained together until 
hydrogenation resulted in the formation of CO2, as sug- 
gested by Goumans & Andersson (2010). In such a case, 
the CO:C02 ratio could be explained by the relative gas- 
phase abundances of CO and atomic oxygen, resulting 
from depletion, while water-ice abundances could remain 
low. However, such a depressed gas-phase oxygen abun- 
dance is likely to be a transient feature of the chemistry 
prior to even greater depletion, as shown by figure 2(b), 
while the severely depressed gas-phase atomic hydrogen 
abundances required to allow oxygen sufficient time to 
diffuse to the nearest CO molecule on the grain surface 
without being hydrogenated are unlikely to pertain in 
the moderate density conditions of dark clouds. 

We therefore suggest that the observed apolar charac- 
teristics of CO2 (and indeed of CO) may not be indicative 
simply of water-poor ice; rather, they would indicate that 
the CO2 is formed from a highly-segregated CO:H 2 
mixture. In such a scenario, CO molecules would form 
islands within a predominantly water-ice composition. 
The formation of OH on these CO-rich islands would 
lead to C0 2 formation on the island, while the immo- 
bility of any OH formed on a water surface would ulti- 
mately lead to its hydrogenation, producing more water. 
This structure would be incorporated into the deeper 
ice layers as new surface layers were formed, produc- 
ing "lumps" of CO within the H20-dominated mantle. 
Oberg et al. (2009c) considered the segregation of ice 
mixtures in protostellar envelopes, caused by the prefer- 
entially higher binding of molecules to each other than to 
different species. They suggested a barrier to CO:H 2 
segregation of 300 ± 100 K, although the relative impor- 
tance of surface or bulk diffusion or swapping processes 
was not well defined. In our dark-cloud regime, the CO 
and H 2 molecules would segregate as the ice mantle 
formed by deposition of material from the gas phase. 
While CO is not sufficiently mobile at around 10 K for re- 
action (2) to occur through direct diffusion, it is mobile 
enough to find other surface CO molecules (which are 
relatively abundant, at around 10 % surface coverage), 
before a new layer of ice is formed and the structure is 
more firmly locked in. Even assuming a dust temperature 
of 9 K, with surface diffusion barriers as high as 0.4 Ed es 
(giving Edif (CO) = 460 K), grain-surface CO molecules 



could (under dark cloud conditions) manage around 10 
thermal hops before a new layer of ice formed - just suffi- 
cient to segregate effectively. As atomic oxygen has to be 
on top of a CO molecule for CO2 to form, all CO2 formed 
at low temperatures is likely to be adjacent to many CO 
molecules, producing the observed spectral characteris- 
tics. Under the high-temperature regime, CO moves to 
OH, giving an essentially random arrangement of CO2, 
which is thus well mixed with H 2 0. Even at tempera- 
tures as high as 20 K, the surface mobility of CO2 is not 
sufficient to re-arrange itself before a new ice layer has 
been deposited. These suggestions are speculative; how- 
ever, the requirement that the accretion rates be small 
compared to rates of surface diffusion could make them 
difficult to test directly in the laboratory. 

The above scenario is also sufficient to explain the gen- 
erally low abundances of polar CO ice observed in dark 
clouds (e.g. Gibb et al. 2000); large islands of CO sur- 
rounded by water ice would produce a predominantly 
apolar signature. 

The collapse models show that, under the simple dy- 
namical model we employ, the first 50 or so layers of ice 
are extremely rich in CO2, while the CO2 formed follow- 
ing the abrupt switch-over to CO dominance is around 
5-8 times less abundant. On this basis, the 15 % po- 
lar CO2 found by Bergin et al. (2005) toward Elias 16 
would suggest a cloud age of around 8 - 9 x 10 5 yr, us- 
ing the riH,max = 4x 10 4 cm~ 3 model. The model CO2 
profiles are also in general agreement with the collated 
ice-observation data presented by Gibb et al. (2000), in 
which polar CO2 is dominant over apolar (although com- 
parison with the other detected ice species suggest an age 
a little over 1 Myr). The collapse models would suggest 
that most of the CO2 observed toward high-extinction 
objects is formed during an early period of low extinction. 
A more extended low-density period would also improve 
the comparison with CH4 ice abundances. The free-fall 
collapse models that we use are naturally very approx- 
imate; more detailed dynamical models would likely re- 
produce specific features of the ice composition more ac- 
curately. 

A further consequence of the C02-formation scenario 
we suggest above is that different grain-size populations 
could exhibit different CO:C02 ratios due to their dif- 
fering temperatures, which should scale as a^ 1 ' 6 (see 
section 3.2). Indeed, due to the sharp switch-over be- 
tween CO2 regimes at around 12 K, one might expect 
a bi-modal distribution of polar/apolar characteristics 
with grain-size, assuming static dark cloud conditions. 
However, the changing dust temperatures engendered by 
the collapse process would complicate this picture, as 
each grain population would acquire both polar and ap- 
olar ices during their evolution. It is also not clear how 
the chemistry would be affected by differing grain sizes, 
as not only the dust temperatures but their accretion 
cross-sections and surface diffusion rates would also be 
altered. The threshold temperature of <~ 12 K could be 
affected by such changes. These questions can only be 
addressed by the use of more comprehensive gas-grain 
chemical models. 

The three-phase model that we employ here allows the 
ice composition formed at any stage of cloud evolution 
to be preserved through later stages. As a result we find 
that, contrary to the findings of Garrod et al. (2006), 
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most carbon ends up as CO, CO2, CH4 or CH3OH, not as 
large hydrocarbons. Because only the surface ice layer is 
active either chemically or physically, the dissociation or 
re-injection into the gas phase of ice material is strongly 
limited, which could otherwise cause hydrocarbon dom- 
inance in the ices. 

The formation of significant quantities of methanol ice 
requires that the deposition of new ice layers be suffi- 
ciently slow for atomic hydrogen to have time to react 
with CO and H 2 CO before they are incorporated into 
the deeper ice mantles. In the dark-cloud case, this oc- 
curs when atomic oxygen begins to deplete substantially 
from the gas phase, reducing the rate of water-ice produc- 
tion, and removing a sink for surface atomic hydrogen. 
This also tends to occur at the same time as CO be- 
comes the dominant gas-phase species other than H 2 , as 
a result of the steady conversion of C and O into CO. In 
the case of a dense collapsing core, we find that very low 
dust temperatures and/or high density conditions lead 
to the rapid deposition of new ice layers; molecular hy- 
drogen deposition becomes especially important in very 
cold, dense conditions. This leaves CO very abundant 
in the outer ice-mantle layers under collapse conditions, 
with high-density cores achieving minimal conversion to 
methanol. Thus, the model suggests either that dense 
and/or collapsing cores that are rich in methanol may 
have spent a significant period at intermediate densities 
and temperatures (at which hydrogenation is more ef- 
ficient), or that the process of star formation has ele- 
vated the dust temperature, such that H 2 coverage is 
insignificant. Such conditions could explain the partic- 
ularly high methanol-ice fractions detected in low-mass 
star- forming cores (Boogert et al. 2008). However, the 
models presented here say nothing about the chemistry 
in such sources where significantly higher temperatures 
or strong UV radiation fields obtain as a result of the 
star-formation process. It should also be noted that the 
activation barriers and tunneling rates for the hydrogena- 
tion of CO and H 2 CO arc somewhat uncertain. 

It is clear that, in order to treat ice chemistry accu- 
rately, models are required that can distinguish the sur- 
face of the ice from the underlying mantle. This is es- 
pecially important for low temperature regimes, where 
mobility within the mantles is likely to be very low, and 
the penetration of surface reactants into the mantle is 
minimal. It is also necessary to use a physically accu- 
rate treatment of the efficiency of barrier-mediated reac- 
tions, and, in the case of C0 2 production, to include an 
indirect formation mechanism to account for the forma- 
tion of OH in the same binding site as a CO molecule. 
Much more accurate determination of various quantum- 
tunneling rates for barrier-mediated reactions would im- 
prove model accuracy. Furthermore, models such as this 
one that can simulate both the grain-surface and gas- 
phase chemistry are necessary to fully understand the 
evolution of dust-grain ices in interstellar environments; 
the accretion rates of gas-phase species strongly affect 
the grain-mantle chemical composition. For an accurate 
treatment of C0 2 ice formation in dark clouds, the inclu- 
sion in the model of some form of collapse process also 
appears to be necessary. 

5. CONCLUSIONS 



We list the main conclusions of this modeling study 
below: 

1. The reaction CO + OH ^ C0 2 + H is shown by 
our models to be efficient enough to account for the 
vast majority of C0 2 formation. If OH is formed 
on top of a CO surface, the only surface diffusion 
required to produce C0 2 is that of atomic hydro- 
gen. 

2. At low temperatures, the grain-surface formation 
of OH on top of a CO molecule allows C0 2 to be 
formed in a ratio CO:C0 2 ~ 2 - 4, under dark cloud 
conditions. These ratios are robust over a wide 
range of temperatures below ^12 K. We identify 
this component with observed apolar CO and C0 2 
signatures in dark clouds. 

3. At higher temperatures (> 12 K), CO is sufficiently 
mobile to compete effectively for reaction with OH, 
producing ice mantles rich in C0 2 and very poor 
in CO. We identify this component with observed 
polar C0 2 signatures in dark clouds. 

4. Visual extinction-dependent chemical models 
appear capable of reproducing the observed 
extinction-threshold behavior of H 2 0, C0 2 , and 
CO. 

5. The decline in dust temperature caused by the ex- 
tinction of the interstellar radiation field at higher 
visual-extinctions causes the change in CO and 
C0 2 chemical behavior on the grains. A sharp 
switch-over occurs in the CO:C0 2 ratio at a tem- 
perature of around 12 K. The visual-extinction 
threshold of CO is thus tied to this temperature 
threshold. 

6. The water-ice extinction threshold is well repro- 
duced with a model gas density of a few x 10 3 cm~ 3 . 
This threshold is associated with the rate of pho- 
todesorption from the grain surface. C0 2 ice for- 
mation occurs in tandem with that of H 2 0, so long 
as CO is present in the gas phase. 

7. The observed polar C0 2 is associated with early- 
time formation in tandem with water ice. We 
suggest that the apolar component is associated 
with C0 2 that exists within a highly segregated 
CO:H 2 ice. This segregation would occur during 
the formation of each ice layer, producing CO-rich 
lumps within the H 2 0. Formation of OH upon the 
CO surface would produce C0 2 that is trapped on 
that surface and would be incorporated into the 
CO-rich ice. A small polar CO signature could be 
produced at the interface of the CO with the H 2 0. 

8. Variations in dust temperature across the canonical 
grain-size distribution may be important in deter- 
mining polar/apolar ice ratios, due to the sharp 
temperature dependence of CO/C0 2 -formation 
rates on the grains. 

9. In the dark-cloud models, large methanol-icc abun- 
dances with respect to water are indicative of the 
depletion of gas-phase oxygen atoms at late times. 
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These large methanol abundances reside in the 
outer layers of the ice mantle. 

10. Under high-density and high-extinction conditions, 
extremely cold temperatures lead to high surface 
coverage of molecular hydrogen, resulting in less 
methanol production due to the shorter times avail- 
able to hydrogenate before new layers are de- 
posited. Thus, ice mantles in cold, dense, collaps- 
ing cores are CO-rich in all but the deepest layers. 
The efficiency of methanol formation at high den- 
sities is extermely sensitive to timescale and dust 
temperature. 

11. Gas-phase methanol is solely produced (on the 
grains) during the later, apolar CO/CO2 phase, 
which occurs after the CO extinction threshold 
has been exceeded. The gas-phase abundance is 
highly stable over this period. 



12. Methane (CH4) and CO2 in the ice mantles appear 
to be well correlated, due to the coincidence of high 
gas-phase atomic carbon abundances (in the case of 
CH 4 ) and elevated dust temperatures (in the case 
of CO2). Methane and CO2 ices are thus mainly 
associated with the deeper ice layers in the mantle. 



The very low observed CH4 ice abundances are in- 
dicative of long periods of conversion of C to CO in 
the gas phase, at visual extinctions in the range 1 - 
2 (from cloud edge to center) , prior to any collapse. 

13. More detailed gas-grain models that deal with the 
short-lived intermediates of CO2 production would 
allow a closer investigation of the relevant reaction 
processes and would improve model accuracy. The 
incorporation of a grain-size distribution into the 
models will also be important to the explanation 
of observational results. 
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Fig. 1. — Growth of ice mantle over time, for static dark-cloud model at T gas = Td us t = 10 K. 
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Fig. 2. — Basic dark-cloud model at 10 K; a) Fractional ice composition; b) Gas-phase abundances with respect to total hydrogen. 
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Fig. 3. — Ice-mantle composition by layer, with dust-grain temperatures: a) 8 K; b) 12 K. Gas temperature is held at 10 K for all runs. 
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Fig. 4. — Ice-mantle composition by layer, assuming -Edif = 0.3E^ CS , with dust-grain temperatures: a) 8 K; b) 12 K. Gas temperature 
is held at 10 K for all runs. 



TABLE 1 
Initial Elemental Abundances 



Species, i 


n(i)/n H a 


H 


5.0(-5) 6 


H 2 


0.5 b 


He 


0.09 


C 


1.4(-4) 


N 


7.5(-5) 


O 


3.2(-4) 


s 


8.0(-8) 


Na 


2.0(-8) 


Mg 


7.0(-9) 


Si 


8.0(-9) 


P 


3.0(-9) 


CI 


4.0(-9) 


Fe 


3.0(-9) 



a A(B) = A B 

b Static dark-cloud models only 



Formation of Interstellar CO2 and Other Ices 



Time (yr) 




10= 



Time (yr) 



(b) 



10 7 




200 



Fig. 5. — Ice-mantle composition by layer, assuming njj = 2 X 10 cm , panels as per figure 2. 
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Fig. 6. — Ice-mantle composition by layer, assuming njj = 2 X 10 5 cm 3 , panels as per figure 2. 




Fig. 7. — Calculated dust temperature through a dark cloud, assuming a = 0.1/im, as a function of edge-to-center visual extinction. 
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Fig. 8. — Ice column densities integrated through a dark cloud of fixed density, assuming a diffusion barrier height E^ ca = 0.3i?dif- j4v 
value represents the edge-to-edge visual extinction, i.e. twice the edge-to-center value used within the chemical code. The three curves 
for each species correspond to timescales of 0.5, 1, and 2 Myr; higher curves represent later times. Each panel shows the model run at a 
different density: a) nn = 2 X 10 3 cm -3 , b) «h = 6 X 10 3 cm -3 , c) nn = 2 X f0 4 cm -3 



TABLE 2 

Binding Energies of Selected Species 



Species, i 


E dcs (i) [K\ 


H 


450 


H 2 


430 


C 


800 


N 


800 





800 


CO 


1150 


HCO 


1600 


OH 


2850 


H 2 


5700 
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Fig. 9. — Ice column densities integrated through a dark cloud of fixed density njj = 6 X 10 3 cm -3 , assuming a diffusion barrier height 
a ) ^des = 0.35-Edifi b) -Edcs = 0.4B(jif; c) E^ es = 0.45i?dif; d) E^ ca = 0.5 i?dif- -Av value represents the edge-to-edge visual extinction, 
i.e. twice the edge-to-center value used within the chemical code. The three curves for each species correspond to timescales of 0.5, 1, and 
2 Myr; higher curves represent later times. 



TABLE 3 

Ice composition, basic dark cloud model" 



Species 


Elias 16 6 


1 X 10 5 yr 


2 X 10 5 yr 


5 X 10 5 yr 


1 X 10 6 yr 


2 x 10 6 yr 


Collapse model, 1.2 X 10 6 yr 


H2O 


100 


100 


100 


100 


100 


100 


100 


CO 


25 


5.8 


10 


16 


20 


20 


30 


CO2 


18 


3.3 


5.7 


9.2 


10 


9 


30 


CH 4 




37 


30 


19 


11 


12 


5.1 


CH3OH 


< 3 


0.34 


1.1 


4.1 


9.7 


18 


8.7 


H 2 CO 




1.2 


2.8 


6.4 


9.3 


10 


12 


NH 3 


< 9 


25 


25 


24 


22 


19 


25 



a Obscrvational values for Elias 16 as collated by Gibb et al. (2000) 
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Fig. 10. — Ice column densities integrated through a dark cloud of fixed density rijj = 6 X 10 3 cm -3 , assuming a diffusion barrier height 
a ) ^dcs = 0.3 -Edit; b) -Edos = 0.4-Edif- Ay value represents the edge-to-edge visual extinction, i.e. twice the edge-to-center value used 
within the chemical code. The three curves for each species correspond to timescales of 0.5, 1, and 2 Myr; higher curves represent later 
times. 
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